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Abstract 

We consider Kaluza-Klein dark matter from latticized universal dimensions. We motivate 
and investigate two different lattice models, where the models differ in the choice of boundary 
conditions. The models reproduce relevant features of the continuum model for Kaluza- 
Klein dark matter. For the model with simple boundary conditions, this is the case even 
for a model with only a few lattice sites. We study the effects of the latticization on the 
differential flux of positrons from Kaluza-Klein dark matter annihilation in the galactic halo. 
We find that for different choices of the compactification radius, the differential positron 
flux rapidly converges to the continuum model results as a function of the number of lattice 
sites. In addition, we consider the prospects for upcoming space-based experiments such as 
PAMELA and AMS-02 to probe the latticization effect. 
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1 Introduction 



Today, cosmological measurements by, e.g., the Wilkinson Microwave Anisotropy Probe 
(WMAP) [1,2], the Two-Degree Field Galaxy Redshift Survey (2dFGRS) [3], and the Sloan 
Digital Sky Survey (SDSS) [4] indicate that the energy and matter in the Universe could be 
distributed such that approximately 73 % is "dark energy" (maybe due to the cosmological 
constant A in Einstein's field equations) and approximately 23 % is "dark matter", which 
leaves around 4 % as ordinary luminous matter in terms of baryons such as protons and 
neutrons, i.e., baryonic matter. 

Further evidence for the presence of dark matter comes from, e.g., galactic rotational curves 
[5,6] and studies of the observed redshift-luminosity relation of type la supernovae [7,8]. 
All these data consistently point towards approximately 23 % dark matter. However, the 
nature of the dark matter is still unknown. Nevertheless, several dark matter candidates 
exist. Ordinary matter could make up a fraction of the dark matter in the form of massive 
compact halo objects (MACHOs) [9,10]. This could be, e.g., faint stars or stellar remnants 
such as black holes. However, the MACHOs can not make up all of the dark matter. 

One of the most plausible candidates for dark matter is weakly interacting massive parti- 
cles (WIMPs) [11, 12]. These are hypothetical non-baryonic particles with masses in the 
GeV-TeV range. They are electrically neutral and assumed to carry a conserved quantum 
number to ensure their stability. Within this category, the prime candidate is the lightest 
supcrsymmctric particle - typically the neutralino [13]. It is electrically neutral and stable 
due to the conservation of R-parity. 

Recently, an interesting alternative dark matter candidate has been intensively studied in the 
hterature, which is so-called Kaluza-Klein dark matter [14,15]. In models with universal 
extra dimensions [16, 17] (for the first model with TeV sized extra dimensions, see [18]), 
in which all standard model particles can move, the first excited state or mode of the 
hypercharge gauge boson is an excellent candidate for dark matter. The prospects for both 
direct and indirect detection are generally good. In particular, for indirect detection, the 
prospects are better than for supersymmetric particles. This is because the Kaluza-Klein 
gauge boson is a vector particle, whereas the neutralino is a Majorana particle, and thus, 
has a hclicity suppressed annihilation cross-section. Therefore, one could obtain an excess 
of positrons from Kaluza-Klein dark matter annihilation in the galactic halo [19,20] (for 
neutralino dark matter detection with positrons, see also [21,22]). This positron flux could, 
e.g., be detected in balloon-borne experiments such as HEAT [23] or future space-based 
experiments such as PAMELA [24] and AMS-02 [25,26]. 

In general, a major disadvantage for models with extra spatial dimensions is that they have 
coupling constants with negative mass dimension and are therefore not renormalizable. 
Hence, these models can only be considered as effective models, valid up to some cutoff 
energy scale, where some more fundamental theory is expected to describe physics. The 
cutoff procedure leads after dimensional reduction to an effective four dimensional descrip- 
tion, but without the full higher-dimensional gauge invariance. One possibility to remedy 
these problems is offered by deconstructed [27] or latticized [28] extra dimensions.^ 

^Note that what is referred to as deconstructed or latticized extra dimensions may differ in the hterature. 
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We will in this paper consider latticized extra dimensions. In these type of models, the 
fundamental high-energy theory is four dimensional with a replicated gauge symmetry of 
the form G x G x . . . x G, where G is a gauge group. In addition, it contains a set of scalar 
fields, so-called link fields, which "link" the gauge groups. If the potential of the scalar fields 
is chosen in an appropriate way, then one can generate the physics of an extra dimension 
in the infra-red regime. This happens when the link fields acquire vacuum expectation 
values and spontaneously break the original symmetry down to a diagonal subgroup of the 
symmetry. The resulting low-energy model is a model with an extra dimension, where 
the extra dimension has been put on a lattice. Note that the model with deconstructed 
dimensions is structurally similar but replaces the link scalar fields with fermion condensates. 

In this way, models with deconstructed or latticized extra dimensions can have similar 
benefits as models with continuum extra dimensions, but since they in contrast with the 
latter preserves manifest gauge invariance and sometimes are renormalizable quantum field 
theories, they have a better defined short distance behavior. This could be relevant to 
models for Kaluza-Klein dark matter, where corrections from physics above the ultra-violet 
cutoff energy scale could be crucial when determining the nature of the lightest Kaluza- 
Klein particle [17]. Radiative corrections to Kaluza-Klein masses have been calculated for 
latticized abelian [29] and non-abelian [30] gauge theories compactified on x S^. So far, 
these calculations have not been applied to the case of latticized universal extra dimensions, 
where there should be additional contributions to Kaluza-Klein masses from the orbifold 
fixed points. In this paper, we will not consider radiative corrections, but instead consider 
the prospects for observing the latticized model already at tree-level. 

More specifically, we will consider two lattice models for a universal extra dimension, where 
the models differ in the choice of boundary conditions. The models are designed to mimic 
in the infra-red regime the continuum theory for Kaluza-Klein dark matter. In the ultra- 
violet regime, the higher-dimensional theory is replaced with completely four dimensional 
dynamics. We study the effects of the latticization and whether one could probe these 
effects in future experiments such as PAMELA and AMS-02. 

The paper is organized as follows. In Sec. |21 we describe the model for a latticized universal 
dimension with simple boundary conditions. Next, in Sec. El we study the flux of positrons 
coming from Kaluza-Klein dark matter annihilations in the galactic halo and the prospects 
for probing the latticization in upcoming space-based experiments. Then, in Sec. |3J we 
consider a lattice model with modified boundary conditions. Finally, in Sec. El we summarize 
our results and present our conclusions. 

2 Latticized Universal Dimensions 

In this section, we will follow closely the results of Refs. [31,32]. 

The model we will consider is a field theory in four dimensions with a product gauge group 
G = n^QS'[/(3)j X SU{2)j X U{l)j. The model contains fermions and gauge bosons as well 
as a set of scalar link fields Qjj+i, ^jj+i, and (pjj+i, where j = 0, 1, . . . , — 1. The link 
fields transform as bifundamentals under adjacent gauge groups. 

Here we follow the classification of Ref. [31]. 
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When the hnk fields acquire vacuum expectation values (VEVs), i.e., {Qjj+i) = "^sls, 
= V2I2, and = Vi/\/2, where 1„ is the nx n identity matrix, the product 

gauge group is spontaneously broken down to the diagonal subgroup SU (3) x SU (2) x 
which we identify as the standard model (SM) gauge group. One can arrange the parameters 
of the scalar potential in such a way that the link fields become non-linear sigma model fields 
in the low-energy effective theory. The low-energy effective theory can then be identified with 
a transverse lattice gauge theory [33], where only the extra dimension has been latticized. 
Since we are interested in a model that mimics the continuum theory for Kaluza-Klein dark 
matter, we will consider a latticized version of an /'L2 orbifold. This topology has in 
the continuum theory the advantages of generating chiral zcroth modes and of removing 
unwanted scalar degrees of freedom^, in agreement with observations. Also, the orbifold 
topology imply a discrete symmetry - Kaluza-Klein parity, which ensures the stability of 
the lightest Kaluza-Klein particle, and thus, making it a viable dark matter candidate. We 
will find that the lattice model considered in this section reproduces the essential features 
of the continuum /'L2 orbifold.^ 

The action for this theory can be split as 5* = 5'gaugc + •S'fcrmion + 5'Higgs- In what follows, wc 
will neglect electroweak symmetry breaking effects and so we will not further consider the 
contribution from the Higgs sector. 



2.1 Gauge Sector 



We will take the action of the gauge sector to be [32] 



i=o 

where a = 1, 2, 3 and b = 1, 2, . . . , 8. Here V{(f), Q) is a suitably chosen scalar potential. 
The Qjj+i fields transform under SU{3)j x S'?7(3)j+i as (3,3), ^jj+i transform under 
SU{2)j X SU{2)jj^i as (2,2), and are charged under U{l)j x [7(1)^+1 as {-Y^,Y^). 

We will set = 1/3. 

The covariant derivatives act on the hnk fields as 



(2) 



= {d, + igA'^^T^^ - i^A^^.+,)^7;'Vi)$,,,+i, (3) 



^The zeroth mode of the fifth component of the higher-dimensional gauge field, ^15°''. 
^We will not explicitly show the absence of A^^\ For a discussion of this, see for example Refs. [28] 
and [32]. 
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where (a = 1, 2, 3) are the generators of SU{2)j, Tj* {b = 1, 2, ... 8) are the generators 
of SU{3)j, and gy, g, and g are the couphng constants of U{l)j, SU{2)j, and SU{3)j, 
respectively. We have assumed discrete translational invariance by setting gyj = gy, gj = g, 
and gj = g, for all j. When the hnk fields acquire VEVs, the kinetic terms for the hnk fields 
generate mass matrices for the gauge bosons. Thus, for the C/(l) gauge fields, we obtain the 
mass terms 



1 2 2 2 ^"' 



(5) 



i=o 



Similarly, for the other gauge fields, we obtain the same kind of mass terms. 
Explicitly, the (A^ + 1) x (A^ + 1) tridiagonal mass-squared matrix reads 



( 1 

-1 



V 



-1 

2 



-1 



-1 



2 

-1 



-1 
1 



(6) 



which can be diagonalized by a change of basis 



N 



(7) 



n=0 



where 



^cos(^7n), n^O 
n = 



1 



Here jn = nn/ (A^ + 1). In this basis, we now obtain 



1 ^ 

n=0 



where 



2 ~2 2-5^2 • 



nn 



(9) 



(10) 



.2(Ar + l) 

In the limit n N, we find a linear Kaluza-Klein spectrum m„ n/R provided that we 
make the identification ngyViY^/[2{N + 1)] = 1/R. It should be noted that when referring 
to R, we refer to the radius of the extra dimension, whereas some other authors with the 
same notation refer to the size of the extra dimension {i.e., ttR for an 5'^/Z2 orbifold). 



2.2 Fermionic Sector 



2.2.1 Fermions 



Next, we include a set of fermions L*; 



and E", for j = 0, 1, . . . , iV and a 



(from now on generational indices will be suppressed). Here Lj transforms as 2 under SU (2)j 
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and as a singlet under SU (2)j, for i ^ j. Furthermore, Lj is charged under U{l)j as Yd = —1. 
The field Ej is a singlet under all SU{2) groups and is charged under U{l)j as Yg = —2. 
Both Lj and Ej transform trivially under all SU{3) groups. 

In Ref. [31], the latticized action for the quark sector was considered. We construct anal- 
ogously the action for leptons as Sfermion = Sd + Sg, where Sd refers to the part containing 
the SU{2) doublet fields and Ss contains the SU{2) singlet fields. Here Sd is given by 



Sd 



N 



V2 [Vi 



+ h.c. 



(11) 



where Mf is a mass parameter which will be used when matching to the continuum model. 
The covariant derivative is given by D^j^Lj = {d^ — igAj^T^ — i^y^A^^) Lj. If the index of 
a field is out of bounds, then it is implicit that the field is zero. For the SU{2) singlet fields, 
we have 



[j=o j=o L ^ 



\{vi/V2r 



E 



E. 



+ h.c. 



, (12) 



where the covariant derivate is given by D^j^Ej = (9^ — igy^Aj^^) Ej. Note that the signs 
for the mass terms in Eqs. ()11|) and (fT^ are different. This is a characteristic feature of 
theories with continuous universal extra dimensions after dimensional reduction [16]. Here, 
as in Ref. [34], it has been obtained by choosing different signs for the Wilson terms when 
deriving Eqs. (fTT|) and p2j] . The signs for the mass terms are also in agreement with the 
corresponding results for the quark sector in Ref. [31]. 

In order to obtain chiral zeroth modes, we take the doublet fields to satisfy Lo_r = and the 
singlet fields to satisfy Eql = 0. When the link fields acquire universal VEVs, we obtain 
mass matrices for the fermion fields. For example, in the doublet sector, we have 



-Mf 



N 



LqlLir + LjL (Lj^i^n — Lju) + h.c. 



(13) 



The corresponding mass matrix is diagonalized by a change of basis 



N 



N 



LjL — ^ ] Q-jnLnL , LjR — ^ ^ hjnLnR , 



(14) 



n=0 



n=l 



where the a,„'s are given in Eq. (jHJ and 



sin Uln) 



N+1 

Thus, we obtain for the left-handed fields the masses 



mi = 4M; sin^ 



rnr 



2{N+1) 



(15) 



(16) 
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where n = 0, 1, . . . , and for the right-handed fields the same form for the masses, but now 
n = 1,2, . . . , N. Thus, there are no zeroth modes for the right-handed doublet fields. For 
n -C A^, we find a linear Kaluza-Klein spectrum if we make the identification 7rMf/{N+l) = 
1/R. Note that with this identification the masses in Eq. become identical to the masses 
in Eq. (jlUp , after the corresponding identification of parameters has been made for Eq. pO|) . 
This is expected and it motivates the use of the same notation in Eqs. (fTUj) and (fTBj) . 

Similarly, we obtain a mass matrix for the singlet sector, which is diagonalized by a change 
of basis 

EjL = bjnEnL , = aj„£'„R. (17) 

n=l n=0 

Note that Ljl has the same expansion as Ejfj, whereas Lj^ has the same expansion as Ejl- 
Furthermore, note that we obtain a right-handed zeroth mode and no left-handed zeroth 
modes for the singlet fields. This is because we have interchanged the role of the left- and 
right-handed fields in this case. 

Note that, as in the continuum model, we obtain negative masses for the singlet fields. 
However, this can be remedied by a change of basis En —'y^En, for n = 1, 2, . . . , [16]. 



2.2.2 Fermion Gauge Boson Couplings 



Now, we consider the couplings between fermions and the U (1) gauge bosons. For the SU (2) 
doublet fields, we obtain from the fermionic kinetic terms Ylf=oLji'y^D^Lj, the couplings 



Qy^ ^j=oLj'y'^AjnLj. In the mass eigenbasis, we have for the left-handed fields 

^ffA = 9Y— O'jnajmO'jlLnLAm^i'y^LiL. (18) 

n,m,l j 

Using the orthogonality relations given in Appendix^ we have 

%A = -gy ^J^^ UlA,,YLol + h.c. + . . . , (19) 

where the dots indicate terms that are not directly relevant to our discussion. In the 
continuum model, the corresponding coupling is given by [16, 35] 

Aont. = \gYYdLfA^^h^{l - 75)^^°^ + h.c. (20) 

There is an additional factor 1 /2 in Eq. (PUI) which is not present in Ref . [35] . This is simply 
due to the fact that we have chosen a different convention for the hypercharge assignment. 
If we identify gy/y/N + 1 = gy, then we observe that we have the same coupling as for the 
continuum case. 

For the singlet fields, we obtain from the kinetic terms ^^^^Eji'-^^D^Ej the couplings 
gy^ "^^=0 Ej'y'^Aj^Ej. In the mass eigenbasis, we have for the right-handed fields 

^ffA = 9Y-^ ^ Q'jnO'jmO'jlEnRAmtil^ Em. (21) 

n,m,l j 
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Again using the orthogonality relations, we have 

^ffA = -9Y^^^^kRA,,rEoR + h.c. + . . . , (22) 

where we have redefined the singlet fields En —>■ —j^En, for n = 1,2, . . . , N, in order to 
obtain positive masses in the singlet sector. 

In the continuum model, the corresponding coupling in the singlet sector is [16, 35] 

Aont. = -^^yy.4'^4^)7'^(l + 75)^^°^ + h.c. (23) 

Thus, with the identification gy/y/N + 1 — gy, we have the same coupling as for the 
continuum model. 

3 Indirect Detection of Kaluza— Klein Dark Matter 

In models of continuous universal extra dimensions, the annihilation of Kaluza-Klein gauge 
bosons can lead to an excess of positrons, which could be observed in balloon-borne and 
space-based experiments. The upcoming space-based experiments PAMELA and AMS-02 
will reach considerably higher energies than the balloon-borne HEAT experiment did. The 
precision will also be largely improved. These improvements are due to the large acceptance 
and long exposure time of the new experiments. We focus in this paper on the PAMELA 
and AMS-02 experiments. 

The space-based (or more specifically satellite-borne) PAMELA experiment is planned to 
be launched in June 2006. It has an acceptance of 20 cm^ sr and an exposure time of three 
years. The PAMELA experiment is optimized to measure the cosmic positron spectrum for 
positron energies up to 270 GeV. 

The space-based AMS-02 experiment will be placed on the International Space Station and 
is planned to be launched in 2008. It has an acceptance of 450 cm^sr and an exposure 
time of three years. The AMS-02 experiment is optimized to probe positron energies up to 
400 GeV. 

In this section, we consider the differential positron fiux from Kaluza-Klein dark matter 
annihilations for the latticized model described in the previous sections. We examine the 
latticization effect and the prospects for upcoming experiments to probe this effect. It should 

be noted that Kaluza-Klein dark matter could also be observed indirectly by observing high- 
energy neutrinos or photons from dark matter annihilations or directly from scattering off 
nucleons [15,36]. We leave the consideration of these cases for future work. 

3.1 Annihilation of Gauge Bosons 

Positrons can be produced from several annihilation channels. They can be produced both 
directly and indirectly through secondary decays. Here we will follow Ref. [15] and only 
consider positrons coming from direct e'^e~ production. Positrons can also be produced 
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f 




Figure 1: Kaluza-Klein gauge boson annihilation into SM fermion pairs. Shown in the figure 
are the t- (left) and u-channel (right) diagrams. 



from cascade decays of muons, taus, and heavy quarks. The lowest-order diagrams for the 
annihilation of the U{1) gauge bosons to particle-antiparticle pairs are given in Fig.Q where 
we will consider / = e, i.e., direct e~^e~ production. Note that the fermionic propagators in 
Fig. ^ can be both left- and right-handed. 

The analysis will be almost identical to that of the continuum case, we only have to replace 
the continuum Kaluza-Klein mode masses with the corresponding masses of the latticized 
model. Thus, from the Feynman rules in Fig. [7|in Appendix ^ we obtain for the t-channel 
annihilation diagrams the amplitudes 



Ml,r = -iglRuYPL,R 



{kp - Pp) Y 
{k — pY ~ ""^1 



(24) 



Here u and v denote four component spinors, and denote polarization vectors, and the 
parameters gi r are given by 



gL,R 



gYYd,s 



gvYd^s 



(25) 



For the w-channel diagrams, one obtains similar expressions. The cross-section for the 
continuum process has been calculated in Ref. [14]. In analogy with their result, we have 



,4 4 , 10 (2m? + s) artanh(/5) - ls(3 
K9l + 9r) 



727rs2/32 



where 



1 



Ara{ 



(26) 



(27) 



mi = 2Mf sin {7r/[2{N + 1)]} is the mass in the propagator, and s = E'^^^ is a Mandelstam 
variable. 



The differential positron flux is given by [37, 38 

d<l>,+ 



Y,^-^vBl+ J de/,(e)G'(e,i?), 



(28) 



where the sum is over all annihilation channels, po is the local dark matter density, and the e"*" 
branching fraction in channel i is denoted by 5*+. The mass m^^ = mi is given by Eq. ()in|l . 
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The Green function G{e, E) describes the propagation of the positrons through the galaxy 
and fi{e) denotes the initial positron energy distribution from channel i. As mentioned 
above, we will only consider the flux of positrons coming from direct e~^e~ production. 
Thus, we will only be concerned with the cross section a^+e-- Since we are concerned with 
non-relativistic particles, we can make a series expansion of ae+e-v for small velocities v. 
Thus, we obtain as in Ref. [15] 



7A 



9y 



288tt{N + l)^mj 



{Yl + • (29) 



Here we have written Eq. (|29|) explicitly in terms of the parameters of our original model. 
Recall that we have made the identification gy/ VN + 1 = gy- Note also that this result 
includes a factor of 1/16, which is not present in Ref. [15]. This is because we have a 
different convention for the hypercharge assignment. 

In calculating the differential positron flux, we have used the DarkSUSY package, see 
Refs. [39,40]. In Figs. |2]|ni we present the differential positron flux as a function of the 
positron energy for an inverse radius of 300 GeV, 450 GeV, 600 GeV, 750 GeV and 900 GeV 
respectively. We present the results for latticized models with N = 1 {i.e., two lattice sites), 
N = 2, and = 3. In addition, we give the continuum model results. We have chosen an 
isothermal sphere dark matter distribution^ with a halo size of z/j = 4 kpc and a local dark 
matter density of po = 0.3 GeVcm"^. The propagation parameters are taken from Ref. [41]. 
In Figs.|21El we also make a comparison with a model for the background differential positron 
flux. In all figures, we have chosen "model C" from Ref. [42]. We have included a factor 1/2 
when calculating the differential positron flux, which comes from the fact that the gauge 
bosons always annihilate in pairs [43]. This factor has not always been accounted for in 
previous works, which is why our continuum results may differ by a factor 1/2 from some 
previous results on the continuum positron spectrum. 

In Ref. [31], the bounds from electroweak precision observables on the mass mi of the first 
Kaluza-Klein mode for latticized and continuous universal extra dimensions are discussed. 
Note that for a continuous extra dimension mi equals the inverse radius R~^, whereas for 
a latticized model we have that mi = 2(A^ + l)(7ri?)-^ sin {7r/[2(A^ + 1)]}. In Ref. [31], 
the most stringent bound for a single continuous universal dimension was found to be 
mi = R-^ > 400 GeV (95 % CL). It was also found in Ref. [31] that the bound for a 
few-site lattice model can be lowered by 10 %-25 %, which would allow for a bound as low 
as mi > 300 GeV (95 % CL) for a few-site lattice model. The reason why the bounds are 
less restrictive for the lattice model is the realization of only a few Kaluza-Klein modes. 
This general feature of lattice models could be important for the PAMELA and AMS-02 
experiments, which for a continuum model could be out of range of the acceptable energies. 
In a recent paper, Ref. [44], the authors argue that the bounds from electroweak precision 
observables could be improved, giving a bound as severe as mi = > 700 GeV (99 % CL). 
This bound was obtained by taking into account two-loop effects and LEP2 data. However, 
as argued above, the bounds for a lattice model should be less severe, especially for a few-site 
lattice model. There are also bounds from WMAP for the relic density of Kaluza-Klein dark 
matter [45]. It was found in Refs. [46,47] that in the minimimal continuum UED model, 

''The Frenk-Navarro- White (FNW) distribution gives very similar results. 
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including coannihilations, to account for the observed relic density, the ideal mass range 
for the lightest Kaluza-Klein particle should be 500 GeV— 600 GeV. Note that the bounds 
from WMAP can be lowered by assuming that Kaluza-Klein dark matter only makes up a 
fraction of the cold dark matter in the Universe. Of course, in this case, the annihilation 
rate in the halo will be reduced, making the detection of positrons with PAMELA and 
AMS-02 more difficult. In addition, there are bounds from direct detection [48], where the 
most stringent bound for a continuum model gives mi > 400 GeV [49], which is no more 
restrictive than the bounds from electroweak precision observables. 

Note the peak in the positron spectrum at a positron energy equal to the WIMP mass {i.e., 
the mass of the first excited mode of the U{1) gauge boson, Ef,+ = m^^ = mi). This peak 
is due to the monoenergetic positron source and is a characteristic signature of Kaluza- 
Klein dark matter. In contrast, for neutralinos, the positron spectrum would be much 
smoother [50]. This is due to the Majorana nature of neutralinos, which leads to helicity 
suppression of direct e~^e~ production. Note that since the lattice model only differ from 
the continuum model in the masses of the Kaluza-Klein modes, one could not by measuring 
the differential positron fiux distinguish between a lattice model with a given radius and a 
continuum model with a larger radius, i.e., both models would make the same prediction.^ 
Such a degeneracy would also exist between a lattice model with parameters A^i and -Ri and 
a different lattice model with parameters N2 and R2. Thus, one would require additional 
input from independent, different types of experiments^ in order to use the information from 
the differential positron flux to probe lattice effects. 

We also observe that the lattice model converges very quickly to the continuum model. 
Given the uncertainties in astrophysical inputs and propagation models, it could be hard 
to detect such small deviations. For increasing values of N, the magnitude {i.e., the height 
of the peak) of the differential positron flux decreases and the peak is shifted towards the 
peak of the continuum model, whereas for decreasing values of the radius R of the extra 
dimension, the separation among the peaks for different values of increases and the 
differential positron flux (as well as the height of the peak) decreases. This can be observed 
in Eq. ()28|) . which shows that the decrease comes from the fact that the differential positron 
flux is suppressed by the mass of the Kaluza-Klein mode. Since the mass of the Kaluza- 
Klein mode is inversely proportional to the radius of the extra dimension, the flux will be 
smaller for a smaller radius, which can also be seen when comparing Figs. I21E1 

4 Modified Boundary Conditions 

Finally, we will consider a model with more complicated boundary conditions. For the gauge 
flelds, we impose that 

- A^„i = 0, (30) 

^At least this conclusion holds at tree-level. We will not investigate radiative corrections. 

^For example, one could compare the level-spacing in the mass spectrum between the zeroth and the 
first Kaluza-Klein modes to the spacing between the first and the second modes. For the continuum model, 
the spacing will be the same, whereas for the lattice model the spacing will decrease. Such an experiment 
would also uniquely determine N and R. 



10 




Figure 2: The differential positron flux (above background) for an inverse radius of 300 GeV as 
a function of positron energy, where we have only considered direct e'^e~ production. In addition, 
positrons can also be produced from cascade decays of muons, taus and heavy quarks. Presented 
are latticized models with two lattice sites (A'' = 1, dotted curve), three lattice sites {N = 2, 
dashed curve) and four lattice sites {N = 3, dash-dotted curve) as well as the continuum model 
(solid curve). Given is also an estimated background flux (gray shaded). The unit of the ordinate 
is cm~^s~^sr~^GeV^-^. 
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Figure 3: The differential positron flux (above background) for an inverse radius of 450 GeV as 
a function of positron energy, where we have only considered direct e+e" production. In addition, 
positrons can also be produced from cascade decays of muons, taus and heavy quarks. Presented 
are latticized models with two lattice sites (A'' = 1, dotted curve), three lattice sites (A'' = 2, 
dashed curve) and four lattice sites {N = 3, dash-dotted curve) as well as the continuum model 
(solid curve) . Given is also an estimated background flux (gray shaded) . The unit of the ordinate 
is cm~^s~^sr~^GeV^-^. For presentation purposes we have chosen the ordinate to start at 0.002. 
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Figure 4: The differential positron flux (above background) for an inverse radius of 600 GeV as 
a function of positron energy, where we have only considered direct e+e" production. In addition, 
positrons can also be produced from cascade decays of muons, taus and heavy quarks. Presented 
are latticized models with two lattice sites (A'' = 1, dotted curve), three lattice sites (A'' = 2, 
dashed curve) and four lattice sites {N = 3, dash-dotted curve) as well as the continuum model 
(solid curve) . Given is also an estimated background flux (gray shaded) . The unit of the ordinate 
is cm~^s~^sr~^GeV^-^. For presentation purposes we have chosen the ordinate to start at 0.002. 
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Figure 5: The differential positron flux (above background) for an inverse radius of 750 GeV as 
a function of positron energy, where we have only considered direct e+e" production. In addition, 
positrons can also be produced from cascade decays of muons, taus and heavy quarks. Presented 
are latticized models with two lattice sites (A'' = 1, dotted curve), three lattice sites (A'' = 2, 
dashed curve) and four lattice sites {N = 3, dash-dotted curve) as well as the continuum model 
(solid curve) . Given is also an estimated background flux (gray shaded) . The unit of the ordinate 
is cm~^s~^sr~^GeV^-^. For presentation purposes we have chosen the ordinate to start at 0.003. 
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Figure 6: The differential positron flux (above background) for an inverse radius of 900 GeV as 
a function of positron energy, where we have only considered direct e+e" production. In addition, 
positrons can also be produced from cascade decays of muons, taus and heavy quarks. Presented 
are latticized models with two lattice sites (A'' = 1, dotted curve), three lattice sites (A'' = 2, 
dashed curve) and four lattice sites {N = 3, dash-dotted curve) as well as the continuum model 
(solid curve) . Given is also an estimated background flux (gray shaded) . The unit of the ordinate 
is cm~^s~^sr~^GeV^-^. For presentation purposes we have chosen the ordinate to start at 0.003. 
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which is a discretized version of the continuum Neumann boundary condition. This will 
lead to the same type of gauge boson mass matrix as in Sec. |2l but now it will be an x 
matrix instead of an {N + 1) x {N + 1) matrix. Thus, the expansion in mass eigenstates 
will now take the form 

N-1 

Aj = ^ ^ ttjnAn, (31) 
n=0 

where the defined as in Eq. (jSj), but with the replacement + 1 N. The 

mass eigenvalues will now be = QyvfY^ sin^ [mT/{2N)], where n = 0, 1,...A^ — 1. 
The correspondence with the continuum mass spectrum is obtained by requiring that 
TT~gYViY^/i2N) = 1/R. 

We impose as in Ref. [32] discretized Neumann boundary conditions in the fermionic sector 
by requiring that Ljl and Ejji satisfy 

Lnl — Ln-i^l = and E^r — -Eat-i^r = 0. (32) 

We take the Ljji and Ej^ fields to satisfy discretized Dirichlet boundary conditions, i.e., 
we impose that Lqji = L^r = and Eql = E^i = 0. In this way, we obtain, as in the 
continuum model, chiral zeroth modes. The corresponding mass matrix is diagonalized by 
a change of basis 

Af-l 7V-1 

LjL = O-jnLnL , Lju = hjnLnR i (33) 

n=0 n=l 

where the a^n's and hjn?> are given in Eqs. (jHI) and (fT^ . respectively, again with the re- 
placement N + 1 ^ N . We obtain for the left- and right-handed handed fields the masses 
m\ = 4:Mj sin^ [nvr/ (2A^)]. For n <ti N, we find a linear Kaluza-Klein spectrum by making 
the identification t^Mj/N = 1/R. 

Similarly, we obtain a mass matrix for the singlet sector, which is diagonalized by a change 
of basis 

N-1 N-1 

EjL = ^ bjnEnL , EjR = ^ ajnEnR. (34) 

n=l n=0 

We obtain as in Sec. 12.2.11 negative masses for the singlet fields, which can be remedied by 
a change of basis En —'y^En, for n = 1, 2, . . . , — 1. 

4.1 Fermion Gauge Boson Couplings 

Analogously to Sec. 12.2.21 "we derive in this section the couplings between fermion and 
gauge bosons for the model with modified boundary conditions. First, we consider the 
SU{2) doublet fields and find for the left-handed fields the couplings 

CffA = ^yy |Aoii,.IoL7''^L + [(Ai + ^/T/N)A,,,LolY'Lil + h.c] 



Af-1 _ >. 

+ + h.c.) + . . . i, 

1=2 ^ 



(35) 
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where again the dots indicate terms that are not directly relevant to our discussion and 



/ATS 



COS (^7i) COS (^7fc) , MO 



Afc = -^amaNk = { ■ (36) 



N 



^ cos (^7i) , k = 



For small or moderate A^, the couplings in Eq. (|35p deviate from the continuum results, see 
Eq. (|2m). The deviation is encoded in the parameters A^. Each of them goes as after 
one factor l/\fN has been absorbed in the redefinition of the coupling constant. Thus, for 
large number of lattice sites, the contribution from the A^. parameters become negligible 
and we reproduce the continuum results. 

Second, in the singlet sector, we find the following couplings 

CffA = ^yy |Aoii^4i?Y^Oi?^ - [(Ai + VW)il;.4K7^^1i?, + h-c] 
N-1 _ >. 

- ^ (^AiA,,EoRrEiR + h.c.) + . . . i, (37) 

1=2 ^ 

where again we have redefined the singlet fields En —>■ —'~f5En, for n = 1,2, . . . , N — 1, in 
order to obtain positive masses in the singlet sector. Precisely as for the doublet sector, the 
couplings in Eq. (j37j) deviate for small from the corresponding continuum couplings, see 
Eq. (j23p . However, in the large N limit, we recover the continuum results. 

Thus, we observe that for a small or moderate number of lattice sites, we obtain for the 
singlet and doublet sectors additional vertices, which are not present in the continuum. The 
Feynman rules for this model are given in Figs. IHlandlHlin Appendix El The new Feynman 
rules will lead to additional diagrams for the annihilation process considered in Sec. 13.11 
Thus, the positron flux will have a more complicated dependence on the latticization than 
in the model with the simpler boundary conditions, c/., Sec.|21 In particular, the degeneracy 
between the lattice model and a continuum model with larger radius is broken. However, 
we also observe that the analogue of Kaluza-Klein parity is explicitly broken. Note that 
Kaluza-Klein parity is necessary to ensure the stability of the lightest Kaluza-Klein mode, 
which is necessary for it to be a viable dark matter candidate. As the number of lattice sites 
is increased, Kaluza-Klein parity is approximately conserved. However, as we have seen, 
the lattice model converges very quickly to the continuum results and so in the region were 
Kaluza-Klein parity is approximately conserved, the deviance may anyway be too small to 
detect. Therefore, at this level, the model considered in this section is mostly of academic 
interest. A deeper analysis would be required to determine the phenomenological relevance 
of this model, but such an analysis is beyond the scope of this paper. 

5 Summary and Conclusions 

We have considered Kaluza-Klein dark matter from latticized universal dimensions. We 
have studied two different models for latticized universal dimensions, where the models 
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differ in the choice of boundary conditions. We have examined to what extent the models 
can reproduce the continuum results for Kaluza-Klein dark matter. For the model with 
simple boundary conditions, we have found that we can reproduce the essential features 
of the continuum model, such as chiral zeroth modes and the relevant couplings, even for 
a few-site lattice model. We have especially examined the effects of the latticization on 
the differential positron flux from Kaluza-Klein dark matter annihilation and the prospects 
for upcoming experiments such as PAMELA and AMS-02 to probe the latticization effect. 
For the model with simple boundary conditions, we have found that the results would be 
equivalent to a continuum model with a larger radius. Thus, in conclusion, the results 
from such an experiment should be used in conjunction with the result from some different 
independent experiment in order to be able to probe lattice effects. We have also pointed 
out that, since the experimental bounds on latticized universal dimensions are less severe 
than for the corresponding continuum model, the prospects for detection for the lattice 
model are better. This could be important for the PAMELA and AMS-02 experiments, 
since these experiments could be out of range to detect the results of the continuum model. 

For the model with modified boundary conditions, we have found a non-trivial dependence 
on the latticization. For a small number of lattice sites, Kaluza-Klein parity is violated, 
which means that, in this range, the model may not be a viable model for dark matter. For 
a large number of lattice sites, Kaluza-Klein parity is approximately conserved. However, 
due to the fast convergence to the continuum results, the prospects for detecting lattice 
effects in this range will be small. 
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A Orthogonality Relations 

In this appendix, we list some useful orthogonality relations for our formalism. In all of 
these relations, we have assumed that ^ n,m ^ N. In the main text, when diagonalizing 
the mass matrices, we consider in Eq. (jH)) the expansion coefficients 



where 7„ = mr/^N + 1). In addition, we have from Eq. ()15j) the expansion coefficients 




(38) 




(39) 
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Note that the expansion coefficients obey the following relations 

N N 



j=0 i=0 



Thus, we have for 71,771,1 ^ and n,m,l <^ N 



N 



j=o V + 1) 

Note that if one of the indices n, m, and / is zero, then we obtain 

^ 1 

Furthermore, if two are zero and one is non-zero, then we obtain zero. Finally, if all three 
of the indices n, m, and / are zero, i.e., {n,m,l) = (0,0,0), then we obtain 

^ 1 

^ajnajmaji = -y==. (43) 

3=0 ^ 

The corresponding relations for the bj^s will not be relevant for the purpose of this paper. 
B Feynman Rules 

In this appendix, we give the relevant Feynman rules for the models considered in the 
text. The Feynman rules are derived by first going to the mass eigenbasis and the simply 
reading off the rules for the vertices. In Fig. [7[ we give the Feynman rules for the model 
considered in Sec. and in Figs. |H1 and El we give the Feynman rules for the model with 
modified boundary conditions considered in Sec. IH Here Pr^l = 1(1=1= 75) are the chirality 
projectors. 
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Figure 7: Feynman rules for the model with simple boundary conditions. 




i-QY^Akl^PL forfc = 2,3,...,A^-l 



Figure 8: Feynman rules for the doublet sector of the model with modified boundary conditions. 
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Figure 9: Feynman rules for the singlet sector of the model with modified boundary conditions. 
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